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ABSTRACT 


Shipboard  antennas  operate  in  the  vicinity  of  the  ship's  super¬ 
structure  which  consists  of  numerous  metallic  scatterers  of  various 
sizes  and  shapes.  The  far -field  patterns  of  shipboard  antennas  are 
perturbed  by  these  scatterers  and  are  therefore  difficult  to  predict. 

In  this  report  a  method  is  discussed  for  approximating  the  far-field 
pattern  in  the  horizontal  plane  of  a  thin  vertical  antenna  (e.g.  ,  a  whip 
antenna)  near  cylindrical  scatterers  of  arbitrary  cross-section  whose 
physical  heights  are  at  least  that  of  the  antenna.  The  method  is  based 
on  approximating  the  antenna  with  ar_  infinite  line  source  and  the 
scatterer  with  an  array  of  thin  infinite  wires. 

Some  representative  calculated  and  experimental  results  are 
presented  which  generally  show  very  good  agreement  between  the 
predicted  and  measured  patterns.  In  addition,  limitations  of  the  method 
are  discussed  and  the  computer  programs  for  calculating  the  far-field 
patterns  are  explained  in  some  detail. 
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FAR -FIELD  PATTERN  PREDICTION 
FOR  SHIPBOARD  ANTENNAS 


Page  1  of  46 


I.  INTRODUCTION 

Radiation  from  shipboard  antennas  is  complicated  by  the  presence 
of  the  many  varied  objects  that  make  up  the  ship's  superstructure.  This 
environment  influences  the  far-field  pattern  of  an  antenna  as  well  as  its 
terminal  impedance  and  can  affect  the  mutual  coupling  between  any  two 
antennas.  This  report  will  be  concer  ted  only  with  the  first  problem, 
namely  the  effect  of  the  superstructure  on  the  far-field  pattern  in  the 
horizontal  plane  of  a  thin  vertical  radiator,  which  is  a  good  approximation 
for  most  MF  and  many  HF  antennas. 

The  technique  used  to  predict  the  perturbed  far-field  pattern  is 
essentially  an  extension  of  Richmond's  [  1  ]  method  for  predicting  the 
scattering  pattern  of  an  arbitrary  array  of  infinitely  long,  thin,  parallel 
wires.  He  has  shown  that  a  solid  cylindrical  conducting  object  can  be 
approximated,  for  the  purpose  of  pattern  calculation,  by  replacing  the 
object  with  thin  parallel  wires  located  at  its  periphery,  if  there  are  at 
least  five  wires  per  wavelength  and  one  wire  defines  each  corner.  For 
a  given  form  of  the  incident  wave  (e.g.  ,  plane  waves,  or  cylindrical 
waves  which  will  be  used  in  this  report)  complex  currents  are  induced 
on  each  of  the  wires.  These  induced  currents  are  then  assumed  to  radiate 
into  free  space  to  form  a  scattering  pattern  which  can  be  calculated  by  the 
far-field  radiation  integral. 

Essentially  the  same  procedure  is  followed  here,  and  a  line 
source  parallel  to  the  wire  array  is  used  to  approximate  a  thin  vertical 
antenna  as  shown  in  Fig.  1.  A  computer  program  then  gives  the  far- 
field  pattern  of  the  line  source  plus  the  .field  scattered  by  the  wires, 
that  is,  the  approximate  pattern  of  the  antenna  in  the  presence  of  the 
scatterer . 

The  question  naturally  arises  whether  or  not  a  technique  which 
employs  infinitely  long  wires  and  an  infinitely  long  line  source  with  uni¬ 
form  current  can  predict  with  acceptable  accuracy  the  far-field  of  a  finite 
antenna  with  a  non-uniform  current  in  the  presence  of  a  conducting  cylinder 
of  finite  height.  The  answer  is  "yes"  within  certain  restrictions.  Justifi¬ 
cation  for  this  positive  answer  was  obtained  experimentally;  the  correspon¬ 
dence  between  the  predicted  and  measured  patterns  was  generally  quite 
good. 
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Fig.  1.  Mathematical  model  for  the  case  of  a  line  source 

radiating  in  the  presence  of  a  cylindrical  structure. 


II.  RESULTS 

A .  General 

The  results  of  this  research  consist  of  ( 1)  a  group  of  computer 
programs  to  calculate  the  patterns  of  the  antenna  (represented  by  a  line 
source)  and  the  cylindrical  scatterers  (represented  by  wire  arrays), 
and  (2)  an  extensive  set  of  measurements  to  establish  the  range  of  valid¬ 
ity  of  the  calculations . 

The  mathematical  formulation  of  the  problem  is  given  in  Appendix  A. 
The  equations  in  this  appendix  are  the  equations  on  which  the  computer 
programs  are  based.  The  computer  programs  for  use  with  a  moderate- 
size  (IBM  1620)  computer  are  given  and  explained  in  Appendix  B.  Three 
passes  through  the  computer  are  required;  the  output  of  the  first  pass  is 
obtained  on  IBM  cards  and  is  used  as  the  input  to  the  second  pass,  and 
the  output  of  this  pass  is  similarly  used  as  input  for  the  third,  which 
yields  the  desired  results.  A  sample  calculation  is  given  in  Appendix  C 
to  illustrate  the  use  of  these  programs.  The  computer  programs  for 
use  with  a  large  (IBM  7094)  computer  are  given  in  Appendix  D.  During 
the  project  it  turned  out  efficient  to  use  the  smaller  computer  for  initial 
investigations  but  the  larger  for  "production"  runs.  The  "production" 
program  is  being  optimized  further  to  allow  the  generation  of  a  handbook 
of  patterns  with  a  reasonable  use  of  computer  time. 

The  experimental  program  was  designed  to  test  the  calculations. 

Two  types  of  approximations  are  involved  in  the  calculations,  and  the 
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experimental  program  was  designed  to  test  these  independently.  One 
approximation  consists  of  replacing  finite  antennas  and  obstacles  with 
infinite  ones  to  simplify  the  computations.  To  test  this  approximation, 
the  models  (described  in  more  detail  below)  consisted  of  antennas  and 
obstacles  between  two  parallel  plates,  of  which  one  could  be  removed. 

With  the  top  plate  in  place,  the  patterns  should  be  the  same  as  for  the 
infinite  structure  used  in  the  calculations,  except  for  finite  groundplane 
effects.  With  the  top  plate  removed,  the  structure  approximates  the 
practical  finite  case.  Comparison  of  the  two  patterns  shows  the  effect 
of  the  approximation.  The  second  approximation  involves  replacing  the 
solid  scatterer  with  an  array  of  wires  and  the  numerical  approximations 
of  the  programs.  The  errors  involved  here  are  estimated  by  the 
differences  between  the  measured  patterns  with  both  plates  in  place 
and  the  calculated  patterns.  Finally,  the  overall  accuracy  can  be 
evaluated  by  comparing  the  computed  patterns  with- the  measured  patterns 
with  the  top  plate  removed. 

B .  Experimental  Model 

The  model  used  for  the  measurements  is  shown  in  Figs  .  2  and  3. 

At  2  Gc  the  structure  is  approximately  5  wavelengths  in  diameter,  exclud¬ 
ing  the  one -wavelength  long  flare.  The  feed  is  off  center,  being  2\  from 
the  nearest  portion  of  the  flare.  The  styrofoam  dielectric  ring  is  one- 
half  wavelength  thick  and  slighU/  less  than  one-half  wavelength  tall.  The 
height  of  the  ring  was  chosen  so  that  the  plate  spacing  would  be  less  than 
X.  /2.  In  order  to  obtain  an  omnidirectional  pattern  with  the  top  plate  in 
place  and  no  obstacle  present  it  was  necessary  to  adjust  the  top  of  the 
exciting  probe  so  that  it  terminated  about _ L_  "  from  the  top  plate.  With 

the  top  plate  on  an  experimental  check  was  possible  of  the  infinite  case 
itself,  whereas  with  the  top  plate  absent  the  more  practical  finite  case 
could  be  measured  and  checked  against  the  patterns  predicted  and  measured 
for  the  infinite  case. 

C .  Comparison  of  Measurements  and  Calculations 

Several  solid  metallic  scattering  objects  were  tested  with  the 
parallel-plate  structure.  These  included  circular  cylinders  with  diameters 
from  \  /4  to  X  and  square  cylinders  whose  sides  were  either  X  / 4or  \  /  2. 
Patterns  were  measured  and  calculated  for  various  distances  between 
source  and  obstacle  using  these  scatterers.  In  all  cases  good  agreement 
between  the  measured  and  calculated  patterns  was  obtained.  Some  repre¬ 
sentative  results  are  discussed  in  the  following  paragraphs. 

Figure  4  shows  a  cylinder  of  square  cross-section  positioned  a  dis¬ 
tance  d  from  a  vertical  source  and  Fig.  5  a,b,c,d  show  the  associated 
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Fig.  2.  Parallel-plate  device  as  it  would  appear  on  an  outdoor  pattern  range. 


Fig.  4.  Square  cylinder  near  a  line  source. 


far-field  patterns  for  various  values  of  d  with  £  =  X.  /4.  The  omnidirectional 
patterns  shown  as  a  reference  for  the  calculated  patterns  correspond  to  the 
antenna  without  scattering  obstacles.  It  can  be  seen  that  in  general  the 
measured  patterns  change  little  when  the  top  plate  is  removed  and  also  that 
the  theoretical  patterns  are  in  good  agreement  with  the  experimental  ones, 
especially  for  the  smaller  values  of  d.  The  reason  for  the  small  lobe  in  the 
forward  scatter  direction  with  the  top-plate  absent  will  be  explained  shortly. 
The  theoretical  pattern  was  determined  by  representing  the  cylinder  by 
5  wires  on  a  side.  Equally  good  results  can  be  obtained  with  as  few  as 
three  wires,  that  is,  a  wire  at  each  corner  and  one  in  between.  The  effect 
of  the  number  of  wires  on  the  calculated  pattern  is  shown  in  Fig.  6. 

Because  of  symmetry  only  that  part  of  the  pattern  between  zero  degrees 
and  180  degrees  is  shown.  Essentially  the  same  pattern  is  obtained  usinn; 
either  5  or  3  wires  per  side.  If  only  4  wires  total  (i.e. ,  one  at  each 
corner)  are  used,  the  pattern  deteriorates  appreciably.  On  the  other  hand 
using  8  wires  on  a  side  (i.e.,  a  density  of  3  2  wires  per  wavelength)  causes 
a  slight  shift  in  the  pattern  opposite  to  that  when  there  were  only  4  wires 
total.  Also,  increasing  the  wire  density  at  the  corners  has  negligible  effect. 
Consequently,  nothing  is  gained  by  using  an  excessive  wire  density  and 
valuable  storage  space  in  the  computer  is  wasted,  a  fact  that  may  be  of 
considerable  importance  in  cases  where  large  and/or  numerous  scatterers 
are  present.  In  addition,  the  number  of  wires  necessary  to  represent  an 
object  may  be  reduced  by  carefully  eliminating  those  wires  on  the  forward 
scatter  side  of  the  scatterer  as  shown  in  Fig.  7.  The  corresponding  far- 
field  pattern  functions  for  the  three  cases  of  Fig.  7  are  shown  in  Fig.  8. 

The  pattern  for  the  case  where  no  wires  are  omitted  is  almost  identically 
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Fig.  5a.  Far-field  patterns  of  a  square 
cylinder,  f  =  X  /4. 

(a)  d  =  X  /4 . 
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Fig.  5b.  Far-field  patterns  of  a  square 
cylinder ,  I  =  X  /4. 

(b)  d  =  \  / 2. 
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Fig.  5c.  Far-field  patterns  of  a  square 
cylinder  ,  f  =  X  /4 . 

(c)  d  =  3X  /4. 
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Fig.  5d.  Far -field  patterns  of  a  square 
cylinder ,  l  =  X  /4 . 

(d)  d  =  X  . 
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Fig.  6  Effect  of  the  number  of  wires  on  the 
calculated  pattern,  d  ~  \  ,  l  -  X  /4. 


Case  8 


7  .  Wires  omitted  on  the  forward-scatter 
side  of  a  square  cylinder. 


Fig.  8.  Patterns  for  the  configurations  in 
Fig.  7. 
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equal  to  pattern  A  while  patterns  B  and  C  show  differences  mostly  in  the 
forward  scatter  direction.  The  reason  some  of  the  wires  may  be  eliminated 
without  seriously  affecting  the  pattern  is  that  the  currents  of  those  wires 
on  the  forward  side  are  quite  small  compared  to  those  on  the  back,  and 
therefore  they  contribute  insignificantly  to  the  far-field  pattern.  When  the 
square  cylinder  is  rotated  45°  about  its  center,  it  was  found  that  almost 
half  of  the  wires  could  be  eliminated  without  seriously  changing  the  pattern. 
This  situation  is  illustrated  in  Fig.  9. 

Several  sets  of  patterns  were  taken  for  the  square  cylinder  X  /  2  on  a 
side  as  shown  in  Fig.  9  except  that  calculations  were  made  with  wires  on 
all  four  sides.  Again  agreement  between  experimental  and  theoretical 
patterns  was  good  as  shown  in  Figs .  10  a,  b,  c,  d.  It  is  interesting  to  note 
that  the  patterns  in  Fig.  10  and  Fig.  5  are  not  too  dissimilar  for  corres¬ 
ponding  values  of  d  even  though  the  size  and  orientation  of  the  scatterer 
are  not  the  same.  This  fact  is  further  supported  by  the  case  of  a  circular 
cylinder  near  a  line  source,  as  illustrated  in  Fig.  11.  Figure  12  shows 
the  predicted  patterns  for  d  equal  to  X  /  2  and  X  for  cylinders  whose 
diameters  are  X  /2  and  X  .  Thus  it  is  seen  that  source-to-scatterer  distance 
is  the  most  sensitive  parameter  and  that  in  general,  for  small  objects, 
the  size  and  shape  of  the  scatterer  are  important  but  to  a  lesser  degree. 

The. experimental  patterns  in  Figs  .  5  and  10  were  recorded  with  the 
height  of  the  source  approximately  equal  to  that  of  the  scatterer;  these 
patterns  exhibit  a  small  minor  lobe  in  the  forward-scatter  direction  that 
is  neither  predicted  by  theory  nor  measured  when  the  top  plate  is  in  place. 


Fig.  9.  Wires  omitted  on  the  forward-scatter 
side  of  a  square  cylinder  rotated  45°. 
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Fig.  10a.  Far-field  patterns  of  a  square 
cylinder  rotated  45° . 

(a)  d  =  \  /4 . 
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Fig.  10b,  Far-field  patterns  oi  a  square 
cylinder  rotated  45°  . 

(b)  d  =  X  /2. 
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Fig.  10c.  Far -field  patterns  of  a  square 
cylinder  rotated  45° . 

(c)  d  =  3\  14. 
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Fig.  lOd.  Far-field  patterns  of  a  square 
cylinder  rotated  45° . 

(d)  d  =  X  . 
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Fig.  11.  Circular  cylinder  near  a  line  source. 


The  existence  of  this  minor  lobe  is  explained  by  Fig.  13.  The  orienta¬ 
tion  of  the  scatterer  is  the  same  as  that  in  Fig.  9  with  ?  =  \  /4,  d  =  X  /4 
and  the  variable  being  the  cylinder  height  h.  The  source  antenna  height 
hs  was  just  under  X  /  2  and  the  experimental  scatterer  was  a  solid  cylinder. 
It  is  apparent  from  Fig.  13  that  the  small  lobe  in  Figs.  5  and  10  is 
caused  by  diffraction  over  the  top  of  the  cylinder  and  that  this  does  not 
take  place  if  the  cylinder  height  is  on  the  order  of  twice  that  of  the  source. 
Also,  for  the  case  where  the  source  is  taller  than  the  scatterer  consider¬ 
able  energy  is  propagated  over  the  top.  Consequently,  the  method  of  far- 
field  pattern  prediction  described  in  this  report  can  not  be  expected  to 
give  very  accurate  results  for  situations  where  the  antenna  is  taller  than 
the  scatterer.  However,  some  insight  into  the  shape  of  the  pattern  in 
back-scatter  direction  may  be  obtained  depending  on  how  much  smaller 
the  scatterer  is  than  the  source  antenna. 

III.  CONCLUSIONS 

It  has  been  shown  that  the  far-field  pattern  in  the  horizontal  plane 
of  a  vertical  antenna  near  cylindrical  scattering  objects  of  finite  height 
can  be  predicted  with  reasonable  accuracy  by  representing  the  scatterer 
by  an  array  of  infinitely-long  parallel  wires  and  the  vertical  antenna  by 
an  infinitely-long  line  source.  The  far-field  predictions  are  valid  so 
long  as  the  scatterer  is  at  least  as  tall  as  the  antenna  itself  and  at  least 
5  wires  per  wavelength  are  used  to  represent  the  object.  The  additional 
restriction  is  made  that  the  number  of  wires  used  must  sufficiently 
define  the  shape  of  the  scatterer.  Also,  in  situations  where  computer 
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Fig.  12a..  Theoretical  far-field  patterns  of 
a  circular  cylinder. 

(a)  d  =  X  / 2,  D  =  X  /4. 
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Fig.  12b.  Theoretical  far-field  patterns  of 
a  circular  cylinder. 

(b)  d  =  X  /2,  D  =  X  /Z. 
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Fig.  12c.  Theoretical  far-field  patterns  of 
a  circular  cylinder. 

(c)  d  =  X.  ,  D=\/4. 
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Fig.  12d.  Theoretical  far-field  patterns  of 
a  circular  cylinder. 

(d)  d  =  \  ,  D  =  X  / 2. 
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Fig.  13a.  Effect  of  relative  heights  of 

source  antenna  and  scatterer. 
Source  antenna  height  =  X  /  2; 
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Fig.  13b.  Effect  of  relative  heights  of 

source  antenna  and  scatterer. 
Source  antenna  height  =\  /Z. 
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storage  is  at  a  premium  it  is  possible  to  omit  those  wires  on  the  forward- 
scatter  side  where  induced  currents  are  relatively  small  and  only  slightly 
affect  the  far-field  pattern. 

This  technique  for  far-field  pattern  prediction  should  find  consider¬ 
able  use  for  the  preliminary  evaluation  of  shipboard  situations  where  a 
vertical  antenna  of  one  or  more  elements  has  in  its  environment  cylindri¬ 
cal  scatterers  with  arbitrary  cross-sections. 

The  program  presented  here  indicates  the  practicality  of  the  prepara¬ 
tion  of  a  handbook  of  patterns  for  a  vertical  antenna  in  the  presence  of 
scatterers  of  various  sizes  and  shapes  which  would  facilitate  the  optimum 
location  of  many  shipboard  antennas. 
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APPENDIX  A 

MATHEMATICAL  DEVELOPMENT  OF  THE  FORMULAS 
FOR  THE  FAR -FIELD  PATTERNS 


The  problem  of  interest  here  is  that  of  determining  the  far-field 
pattern  in  the  horizontal  plane  due  to  a  vertical  antenna  radiating  near 
conducting  cylinders  of  arbitrary  cross -section.  To  obtain  a  mathematical 
model  we  approximate  the  conducting  cylinder  of  arbitrary  cross-section 
by  an  array  of  infinitely-long  thin  wires  along  the  periphery  and  represent 
the  vertical  antenna  by  a  vertical  line  source,  as  shown  in  Fig.  1. 

The  line  source,  whose  current  is  directed  in  the  positive  z-direction, 
produces  only  a  z-component  of  the  electric  field  which  is  given  by 


(1)  E2A  =  IH0W(kp) 

where 

(2)  k2  =  oj2  (j.€ 

and  HQ(2)  (kp  )  represents  the  zero-order  Hankel  function  of  the  second  kind. 
The  current  induced  on  wire  n  is  uniform  and  produces  a  scattered  field  of 
the  form 


(3)  EnS  =  -  (^)  ln  Jo  (kan)  H0(2>(kPn) 

for  Pn^an, 


where  JQ  (kan)  is  the  zero  order  Bessel  function,  an  is  the  radius  of  wire  n, 
and  pn  is  the  distance  from  wire  n.  For  convenience  a  "modified  current" 
In  is  defined  by 

(4)  In  =  In  Jq  (^an)  • 
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Thus  the  field  scattered  from  N  wires  at  the  surface  of  wire  m  can  be 
written  as 


(5)  EZS  =  -  ^  Xn'  H0(2)(kpmn> 

n  =1 


where  pmn  is  the  distance  between  wires  m  and  n, 


<6>  Pmn  =  ,J(Xm  -Xn)*  +(Ym  -Yn)2 


and 


(7)  Pmm  ~  am  , 


To  obtain  values  for  the  complex  current  on  each  wire,  we  impose 
the  condition  that 

(8)  E^+E^O 

at  the  surface  of  each  wire.  Actually  we  require  that  tangential  E  vanish 
at  the  center  of  each  wire,  which  is  an  acceptable  approximation  for  wires 
whose  radii  are  much  smaller  than  the  wavelength.  Substituting  Eq.(l)  and 

Eq.(5)  with  —  I  =  -  1  into  Eq.(8)  we  obtain 
4 

W>  2  =Hol!)(kPm)  • 

n=l 

The  above  equation  then  yields  a  set  of  N  equations  in  the  N  unknowns  ln'  . 
Having  determined  the  complex  currents  on  each  wire,  the  far-field 
pattern  can  then  be  determined  from 


(10)  Ez<4>)  =  +  ^  H0<2>(kpn)  In*  -  H0(*)(kp)  . 

n  =1 

When  the  observation  point  is  at  a  great  distance  from  the  wires,  it 
is  said  to  be  in  the  far-zone  field.  This  permits  some  simplification  of 
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£q.(lO)  .  For  large  arguments  the  Hankel  function  can  be  approximated 
by  the  asymptotic  formula 


(11)  Ho<z)(X)  =  \|  2j/irX  e‘jX,  X»l. 


The  distance  from  wire  n  at  (Xn,  Y  )  to  the  observation  point  at  (p  ,d>)  can 
be  expressed  as 


(12)  p  n  =  p  -  Xn  cos  <j>  -  Yn  sin  * 


Consequently  the  far-field  is  given  by 

(13) 


Ez  <4>)  =  + 


4 


N 

2i  -jkp 

,  \T  '  j'k(Xncos4>  +  Ynsin4>) 

TT  k  p 

-it  ^in  e 

n  =1  _ 

A  far-field  pattern  function  is  defined  by  the  bracketed  expression  in  Eq.(13), 
which  can  be  written  as 


N 

(14)  F  (4>)  =  ^!n'  ejk(Xncos4»+Ynsin4») 

n=o 

where  I0'  -1  at  the  origin. 


Appendix  B  gives  three  computer  programs  for  implementing  these 
calculations  on  a  moderate-size  (IBM-1620)  computer.  The  first  calculates 
the  coefficients  and  excitation,  functions  for  the  system  of  Eq.(9).  The 
second  solves  the  system  for  the  modified  currents  In’ .  The  last  calculates 
the  far  field  from  these  currents,  according  to  Eq.(14).  If  only  the 
scattered  portion  of  the  field  is  desired,  the  value  of  the  current  on  the 
line  source  is  set  equal  to  zero  in  the  input  to  this  last  program. 

Appendix  C  gives  a  sample  problem  that  illustrates  the  use  of  the 
three  FORTRAN  programs  on  the  IBM  1620. 

Appendix  D  gives  a  single,  more  efficient  SCATRAN  program  for 
use  on  the  IBM  7094  computer,  which  may  be  used  in  place  of  the  three 
FORTRAN  programs. 
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APPENDIX  B 


Figure  14  is  a  computer  program  in  FORTRAN  (Ohio  State  Univer¬ 
sity  Version  2)  based  on  Eq.(9)  for  setting  up  the  linear  equations  for  the 
currents  on  the  infinitely-long  parallel  wires.  The  line  source  is  assumed 
to  have  a  negative  modified  unit  current.  Included  in  the  program  is  a 
subroutine  for  calculating  the  Bessel  and  Neumann  functions.  If  the 
argument,  kpmn,  is  less  than  4  the  ascending  power  series  is  used  a"d  if 
it  is  greater  than  4  the  asymptotic  series  is  used.  The  numbers  in  the 
"equals  column"  represent  the  incident  electric  field  intensity  at  wire  m 
due  to  the  line  source.  Figure  15  gives  a  flow  diagram  for  the  program 
in  Fig.  14.  The  output  of  the  program  in  Fig.  14  is  used  as  part  of  the 
input  data  for  the  program  in  Fig.  16.  This  program  uses  the  method  of 
Crout[2,3]  to  obtain  the  solution  for  an  inhomogeneous  system  of  linear 
equations  with  complex  coefficients.  The  output  of  this  program  consists 
of  the  real  and  imaginary  parts  of  the  current  on  each  wire. 

Figure  17  gives  a  flow  diagram  for  the  program  in  Fig.  16.  However, 
this  program  does  not  compute  the  elements  in  Grout's  "auxiliary  matrix" 
in  the  order  called  for  by  him  in  his  article[2]  .  The  order  used  here 
results  in  a  more  compact  and  efficient  program. 

The  output  of  the  program  in  Fig.  1 6 ,  the  various  wire  locations 
as  used  for  input  to  the  program  of  Fig.  14,  and  the  magnitude  of  the 
line  source  current  are  used  as  the  input  for  the  program  in  Fig.  18. 

This  latter  program  calculates  the  far-field  pattern  function  of  the  parallel 
wire  array  excited  by  a  line  source  whose  "modified  current"  equals 
-  1  +  jO.  If  only  the  scattered  field  is  desired,  the  input  card  for  the 
magnitude  of  the  line  source  current  is  omitted  from  the  input  deck.  This 
program  is  based  on  Eq.(14).  The  output  data  gives  the  value  of  the  pattern 
function  at  whatever  values  of  <j>  are  desired.  Note  that  E(<(>)  in  this  pro¬ 
gram  corresponds  to  the  F(<{>)  defined  by  Eq,(14}. 

Figure  19  gives  a  flow  diagram  for  the  program  in  Fig.  18. 


30 


I 


c 

c 


G  THIELE.  LINEAR  EUUATIONS  FOR  SCATTERING  PATThRN  OF  AN 
ARBITRARY  ARRAY  OF  THIN  //IRES*  LINE  SOURCE.  AT  THE  ORIGIN 
D I  MENS  ION  u < 56 ) *  V ( 56 ) *  F ( 29 ) . G ( 29  ) 

2  FORMAT  (F1J.6*F10.6*F10.S«F1J.6.F10.6.F1J.6*F1J.6.F10.6) 

3  FORMAT  (  I  1 Ui F 1 J«6«F 1  j*6iF  1  O  •  6  «  F 1  L  •  6 «  F  10.6.F10.6.F10.6) 
READ  3tNiBKJt BK Y  .  A  A  A  .  T  HE  T  A 


READ2.U (  1  )  i U ( 2  I  «  U  (  3  )  •  U  <  4  )  «U(5).U(6)  •  U  (  7  )  »  U  l  8  ) 

RE  AD2  i  U  (  9  I  .U(  U ) «U(  1 1  ) .U( 1 2 ) .U( 13 ) .U  < 14 ) .U( 15) .U ( ’ 6 ) 
READ2  >U (  1  7  )  .  U  (  18).UI  19). U< 20). 0(21  )  .U  <22  )  *0(23) 
READ2. V (l).V(2).V(3).V(4).Y(5).V(6).V(7).V(6j 
REA  02  .V(9).V(lO).V<ll).V(12).V(13).V(14).V(l  3*!  .  V  (  1  6  > 
READ2. V ( 1 7  )  . V (  lb  )  . V (  19 ) . V(20  )  . V(2 1  )  . V  < 22  )  » V ( 23  > 

M  =  N 
NN  =  M+l 

THETA  =  .0 1 745329*THETA 

STHET  =  S 1 N ( THE  T  A )  — 

IF (AAA)  8. 5, 3  y 

5  STHET  =  6.2831 853* STHET 
B  DO  10  I.  =  l.M 


1  0 


1  1 


1  b 


1  7 


24 

25 
12 


1  4 


20 


30 


Symbols  for  Input  Data 


?!'J 
B  KY 
AAA 


Number  of 
sin 


wires 


THETA 

U(N) 

V(K) 


l  .M 

1  .  M 


1  .N 


U (  I  )  =  STHET*U  <  I  ) 

V  <  I  )  =  STHET *V <  I  ) 

STHET  s  SlN(THETA) 

DO  lb  I  =  1.24 

F  <  I  )  =  0  . 

G  (  I  )  =  0. 

DO  50  I  = 

00  17  K  = 

F  (  K  )  =  0  . 

G ( <  )  =  0. 

DO  30  J  = 

L  =  J 

IF (M-J 124.2b.25 
L  =  J-M 

IF ( I-J  )  14, 12. 14 
F ( J  )  =  BK J 

G  (  J  )  =  -BKY 
GO  TO  30 

R  =  SQRT (  (U(  I  )-J( J  )  )*(U (  I  ) 

UK  =  1 
GO  TO  34 
F(L)=F(L)+B 
G(L)  =  G(L)-V 
CONT I NUE 

R  =  SORT (U( I >*U( I >+V< I )*V( I > > 
I  JK  =  2 
GO  TO  34 


J0(ka 
H0  (ka 
0 °  if 


°o\ 

°0) 

wire  coordinates 


sin 


1.  if  wire  coordinates 

90. 

x  coordinate  of  wire  N 
y  coordinate  of  wire  N 


are 

are 


x/A  and  y/A 
kx  and  ky 


u  (  J  )  )  +  ( v (  I  ) -v ( j  >  >  *  <  v (  1  ) —v ( J  )  )  ) 


Fig.  14.  Digital  computer  program  for  setting  up  the  linear 
equations  for  the  N  unknown  currents . 
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45  F i NN )  =  B 
G(NN)  -  —  Y 

PUNCH2  .F(l).F(2>.F(3>.F(4>.F(5>.F<6).F(7).F(d> 

PUNCH2  .F(9).F(lu),F<ll).F<12).F(13).F<14>.F(l:j).F(16) 
PUNCH2  <F  (  l  7)  »F  (  18)  .F  (  19)  .F  (20  >  .F  <2.1  >  «F<22  >  .F  (23)  «F  <  24  > 
PJNCH2  <G (l).Gf2).G(3)»G(4)«G(5).G(6)tG(7)«G(Q) 

PUNCH2 . G (9  )  .G(  1 u  )  «G ( 1 1  ) *  G ( 1 2 ) «G  < 13 > .G  < 14  > *G< 15 ) *G(  16 ) 
50  PUNCH2  <G ( 17)»G(18)«G<  19>.G(20).G<21  ).G(22>.G<23).G<24> 
STOP 

34  IF  (  R—  4  •  )  35*35. 40' 

35  W  =  R/2. 

W  1  =i A/*W 

«/2  =  W  1  *wl. 

W3=W2*W1 
W4=W2*W2 
W/5  =  W4*W  1 
W6  =  W3**/3 
\ti  7  =  W  3  *  W  4 


B  =  1  .-Wl  +  .25*rf/2-.027777d**/3+.00  1  7361  *l»/4 
B=3-.6944E-4*W8+. 1 929E-5*W6-. 394E-7*W7 
S  =  W  1-.  375  *v«/2+.  0509259**/ 3-.  00361  69*W4 
S  =  S+.  15d56E-3**.b-.4726E-5*W6  +  .  1021E-6*W7 
Y  =  .63661 97  7* ( S+B* ( . 5772 1 566+LOG ( W  )  ) > 

GO  TO  44 
40  Z  -  l./R 
X-Z*Z 

P= (  1 .-.G7uJl *X+« 1  121 5*X*X  >*.564 1 896 

Q=( . 1 25- • o  7324 *X ) *Z* • 564 1 896 

B= ( (P-U )*COS<R)+(P+Q)*S IN(R ) ) /SORT (R ) 

Y=( (P-0  3IN(R)-(P+0)*COS(R) ) /SORT ( R ) 

44  GO  TO  <  20.45  )  «  I JK 
END 


Fig.  14.  Digital  computer  program  for  setting  up  the  linear 
equations  for  the  N  unknown  currents. 
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Flow  graph  for  the  program  in  Fig.  14. 


C  J  RICHMOND*  SOLUTION  OF  15  (Ok  LtSd)  COMPLEX  LINEAR  EQUATIONS 

DIMENSION  A(  1  b  «  1  6  )  *  d  (  15«16> 

2  FORMAT  (F10.6»Fl0.6«F10.6*F10.6«F10#6*F10.6»F10*6»Fi0.6) 

24  FORMAT  < I5.F15.6.F15.6) 

READ  l.N 
DO  5  I =1 .N 

READ2*  A (l*l).A(I*2)*A(I«3>*A(  1*4)*A(I«5)«A(  1«6)*A(1*7)«A(1*Q) 


READ2  •A(l*9)«A(I*lu)*A(l*ll)*A(l«12)*A(I*13)*A(I*14)«A(l*15)*A(I*16) 
READ2  *  d  (  1  *  l  )  iB(  l  *2  )  *d  (  I  *  3  )  *d  (  1 . 4  >  *0  <  I  *5  )  *d  (  1  *6  )  *B  (  I  «  7  )  *B  <  I  *8  ) 


17 

4 


READ2  *  D  (  1  *  9  )*  13  !  I  *  1  0  )  *  6<  1*11  )  « d  <  I  ♦  1  2  )  *B(  I  «  1  3  >  ♦  8  1  I  *  14  )  »  B  1  l  *1  5  )  «  B  (  1  «  1  6  ) 
NN  =  N  +  l 
DO  18  L  =  1.N 
LLL  =  L-l 
DO  IB  I  =  L«N 
11=1+1 
IF ( LLL >5.4*3 
DO  17  K  =  1 « LLL 
A  <  I  .  L  )  =  A  (  I  .  L  )  - 
8(1«L)  =  l3  (  I  i  L  )  - 

A  1  L  .  I  I  )  =  A  <  L  »  1  1  > 

6<L*  I  1  )  =  U ( L  *  I  1  ) 

R  =  A(L*L>+A(L*L) 

S  =  A (L*  I  1  > 

T  =  6  ( L  *  I  l  >  ■ 


A  (  I  *  K  )  *  A  (  K  *  L  > 

A  (  1  *K  )*d (<*L > 

-  A ( L • <  )*A(K.  II)  + 

-  A<  L.K )*B1K ♦  I  I  >  - 


+  d (  1  .< )*8 (<*L  ) 

-  A  (K  *L ) *8  (  1  *  K  ) 

a< l.k  >*a (k«  i  i  > 

13  (  L.K  >*A  (K.  I  1  > 


+  d<L*L )*d (L *L > 


A  <  L  *  I  I  )  =  (  A  (  L  *  L  )  *3  +  B<L*L)*T)/R 
18  B ( L  *  I  1  )  =  (A(L.L>*r  -  B<L*L>*S)/R 
PUNCH  24  * N *  A  <  N  * NN ) .BIN* NN ) 


DO  23  L  =  2 « N 
I  =  NN-L 
11=1+1 
DO  22  K  =  1  I  *N 

All. NN )  =  A(I.NN)  -  A(  1  *K)*A(K.NN>  +  d (  I  • K  )  *B  CK »  NN ) 

22  Bll.NN)  =  B(l.NN)  -  A  1  l ♦ < > *B < K. NN )  -  d 1  1  . K  ) *A 1 K . NN > 

23  PUNCH  24  *  1  *  A  <  I «  NN ) . d 1  1  «  NN ) 

STOP 

END 


Symbols  for  Input  Data 


Symbols  for  Output  Data 


N  lumber  of  Linear  equations 

Real  part  of  coefficient  of  i'. 
in  equation  i  ^ 

B(I,J)  Imaginary  part  of  coefficient 
of  I j  in  equation  i 


I  Index  of  current  1^ 

A(l,NN)  Real  part  of  current  1^ 

B(I,NN)  Imaginary  part  of 
current  1^ 


Fig.  16.  Digital  computer  program  for  solving  the  system  of 
of  linear  equations  for  the  N  unknown  currents. 
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C  J  RICHMOND*  PATTERN  OF  ARRAY  OF  LINE  SOURCES 
01  MENS  ION  CR ( 30 ) »CI  (3J )  *X (3j ) *  Y (30 ) 

2  FORMAT  ( J 5.F 1 b .b.F 1 5.b. F 1 b.b .F15.8.F1 5.8 > 

3  FORMAT  <Flo»6*Fl'^.6*F10.b*Flij.b*FlJ#6«Fl0.b*Fl>j.6*Fl0«6) 
READ  2 • N. P* DP « FM , AAA . THE T A 

READ  3*X(1).X(2)*X(3)*X(4)*X(5)«X(6)«X(7)*X(B) 

READ  3»X(y).X(l^).X(ll).X(l2)*X<13)*X<l4)*X(lb).X(16) 
READ3*  X ( j  7 ) . X  t  lb  > . X ( 1 y )  *x (2 J  )  .  X  <2 1  )  « X  < 22  )  • X ( 23 ) ♦ X ( 24 ) 
READ  3  <  Y (  1  )*Y(2)*Y(3).Y(4)*Y(5)*Y(6)*Y(7).Y(b) 


READ  3«Y(9)*Y(  1 o  j  . Y  <  1  1  )  «  Y  C 1 2  ) 
RE AD3 «  Y ( 1 7 ) • Y (  18 ) . Y (  1 V ) * Y (20  ) 
THETA  =  .0 1 74S32S*TH£TA 

STHET  =  SIN(THETA) 

I F  <  A AA >8.5.8 

5  STHET  =  6.28318b3*STHET 
8  DO  1 0  I  =  1  .  N 

X (  I  )  =  STHET*X  (  I  ) 

1 0  Y  <  I  )  =  STHET*Y  (  I  ) 

12  DO  20  1  =  1  ,N 

READ  2.J.U.V 
CR  <  J  )  =  U 
Cl  ( J  )  =  V 

cs  =  u*u  +  v*v 

CM  =  SQRT(CS) 

20  PUNCH  2. J.U.V.CS.CM 
M  =  FM 

DO  4  0  I  =  1  ,M 

PR  =  .01745329*P 
C  =  COS (PR) 

S  =  SIN(PR) 

ER  =  0. 

El  =  0. 

DO  30  J  =  1  ,N 

2  =  X(J)*C  +  Y(J)*S 
CZ  =  COS(Z) 

SZ  =  S IN (Z ) 

ER  =  ER  +  CR ( J ) *CZ  -  Cl  <  J  )*SZ 
30  El  =  El  +  CI(J)*CZ  +  CR  <  J  >  *  SZ 
ES  =  ER*ER  +  E I  *E  I 
E  =  SORT (ES  ) 

PUNCH  2. I .P.ER.E I *ES*E 
40  P  =  P  +  DP 
STOP 
END 


Y  (  1  3  )  ♦ 

Y(14)*Y(lb)*Y(16) 

Y  (  2  1  )  . 

Y ( 22 )  *  Y (23 ) *  Y ( 24  ) 

Symbols  far  Input  Data 

N 

Number  of  wires 

P 

Initial  value  of  p  in  degrees 

DP 

Increment  in  0  in  degrees 

FM 

Number  of  values  of  0 

AAA 

0.  if  coordinates  are  x/}».  and  y/). 

1.  if  coordinates  are  kx  and  ky 

THETA 

Angle  of  incidence  0Q  with 

respect  to  z  axis 

X(N) 

x/\  or  kx  (coordinate  of  wire  n) 

TOO 

yfi.  csr  ky  (coordinate  of  wire  n) 

J 

Index  number  of  wire  j 

u 

Real  rr-rt  of  current  on  wire  j 

V 

Imaginary  current  on  wire  j 

Symbols  for  Output  Data 

P 

Angle  0  in  degrees 

SR 

Real  part  of  5(0) 

31 

Imaginary  part  of  E(0) 

ES 

Magnitude  scuared  of  8(0) 

E 

Magnitude  of  E(0) 

Fig.  18.  Digital  computer  program  for  the  far-field  pattern  of 
an  arbitrary  array  of  parallel  wires . 
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APPENDIX  C 


The  purpose  of  this  appendix  is  to  minimize  the  time  required  to 
gain  proficiency  in  the  utilization  of  the  programs  in  Appendix  B.  As 
a  sample  calculation,  the  situation  in  Fig.  4  was  used  with  d  =  X.  /2  and 
£  =  X.  / 4.  Three  wires  were  used  on  a  side  as  opposed  to  the  five  that 
were  used  to  obtain  the  theoretical  pattern  in  Fig.  5b. 

The  input  cards  to  the  program  of  Fig.  14  are  shown  in  Fig.  20 
and  are  self-explanatory  except  possibly  for  the  four  "zero"  cards 
which  are  necessary  because  the  program  can  accommodate  up  to  24 
wires  and  only  8  are  used  in  this  instance.  The  output  from  this  pro¬ 
gram  consists  of  the  real  and  imaginary  parts  of  the  coefficients  of 
^n  .  Again,  since  the  program  provides  for  up  to  24  wires  but  only  8 
are  used  in  this  example,  some  of  the  output  cards  contain  only  zeros. 
These  "zero  cards"  are  removed  before  using  the  output  as  the  input 
for  the  program  of  Fig.  16.  The  output  of  that  computation  is  just 
the  real  and  imaginary,  parts  of  the  currents  on  each  of  the  8  wires. 
These  cards  are  then  used  as  part  of  the  input  to  the  program  of  Fig.  18. 
Note  that  the  same  coordinate  information  used  for  Fig.  14  is  also  used 
for  Fig.  18.  Also  note  that  the  "current"  -  1  +  jO  of  the  line  source 
must  be  included  in  the  input  unless  just  the  scattered  field  is  desired. 
The  output  of  the  Fig.  18  program  contains  the  desired  far-field  pattern 
information  in  the  right-most  column.  These  numbers  may  be  used  to 
generate  the  normalized  pattern  in  Fig.  5b  which  was  calculated  using 
5  wires  per  side  and  therefore  has  slightly  deeper  nulls .  If  in  this 
example  the  only  wire  were  the  line  source  itself,  the  output  of  the  last 
program  (Fig.  18)  would  indicate  that  F($)  =  1.  Hence  the  unit  circle 
may  be  used  as  a  reference. 
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c 

1  NPuT 

8 

CARDS  FOR 
. 9V937 

PROGRAM  IN  FI 
-1.97931 

lo  1  it 

0.0 

90.0 

0.500 

0.0 

0.0 

0.500 

0.625 

0.750 

0.750 

0 • 75  J  ' 

0  •  o d'o 

0.500 

0.000 

0.0 

c.o 

0.125 

0.  125 

0.125 

0.000 

-0.12b 

-0# 125 

-0. 125 

C  OUTPUT 

cards  from  program 

IN  FIG  14 

+.999370 

+.851631 

+.714556 

+.365432 

+.472001 

+.365432 

+.714556 

+.851631 

-.304243 

+.000000 

+  a  0  U  0  U  0 

+ . U 00000 

+.000000 

+. OOOOOO 

+ 1  OOOjuO 

+  .OOOOOO 

+.000000 

+ .OOOOOO 

+ • OuOOJu 

+.000000 

+. OOOOOO 

+. OOOOOO 

+, OOOOOO 

+ .OOOOOO 

+1 .979310 

+. 1 01207 

1 6960 1 

-.467056 

-.4 1 0003 

-.467056 

1 6960 1 

+  .  101207 

-.328366 

+.000000 

+  •  0  \j  0  0  0  0 

+.000000 

+.000000 

+. OOOOOO 

+. OOOOOO 

+  .OOOOOO 

+.000000 

+.000000 

♦•OuOOOO 

+.000000 

+.000000 

+  . OOOOOO' 

+  .0000t'0 

+.000000 

+.851631 

+ .999370 

+.851631 

+.472001 

+.365432 

+.098474 

+.365432 

+.472001 

-.329898 

+ .oooooo 

♦•OuOOOO 

+. OOOOOO 

+. OOOOCO 

+.000000 

+. OOOOCO 

+. OOOOCO 

+.000000 

+.000000 

+ .OUOOOO 

+. O00000 

+. OOOOOO 

+ .OOOOOO 

+.OOOCCC 

+ . oooooo 

+. 10 1207 

+  1  .97-+310 

+• 1 01 2J7 

-.41 0003 

-.467056 

-.520633 

— . 4q7056 

-.410003 

-.292726 

+ .OOOOOO 

+  0  u  u  0  o  0 

+.000000 

+. OUOOOO 

+ .OCCCOO 

+.CCCOUO 

+ .oooooo 

+ .000000 

+. OOOOOO 

♦•OGOOOO 

+.000000 

+.000000 

+. OOOOOO 

+ .OOOOOO 

+.000000 

+.71 4556 

+.651031 

+ .999370 

+.851631 

+.714556 

+.365432 

+.472001 

+.365432 

-.396862 

+.000000 

+  . 3  uuCUO 

+.000000 

+.000000 

+ . OOOOOO 

+.OOOuOO 

+. OUOOOO 

+ . OOOoOO 

+ . OOuOOG 

+.000000 

+.000000 

+• OOOOuO 

+.000000 

+. UOOOOO 

+ • OOOOOO 

-. 169601 

+. 101207 

+1 .979310 

+. 101207 

-. 169601 

- .467056 

-.410003 

- . 467056 

+.018899 

+.000000 

+.000000 

+. OOOOOO 

+.000000 

+. OOOOOO 

+.000000 

+.000000 

+. OoOOOO 

+. OOOOOO 

+. OOOOOO 

+.000000 

+.000000 

+. OOOOOO 

+.000000 

+.000000 

+.365432 

+.472001 

+.851631 

+.999370 

+.851631 

+.472001 

+.365432 

+.098474 

-.2471 16 

+.000000 

+ . OOOOOO 

+. OOOOOO 

+  . OOOOOO 

+ .OOOOOO 

+. OOOOOO 

+. OOOOOO 

.+  •  000000 

+ .OOoOOO 

+ • OGuOOO 

+.DOOCCU 

+ . uuGOOO 

+ .UOUOOO 

+. OOOOOO 

+. OOOOOO 

-.467056 

-.41 0003 

+. 1 Ol 2u7 

+1 .979310 

+. 1 0 1207 

-.4 1 0003 

-.467056 

-.520633 

+.267421 

+ .OOOOOO 

+  .  0  *'0  0  J 0 

+. JJOOOu 

+ . uGOOUO 

+ .UOOOOO 

+. OOOOUO 

+.000000 

+.000000 

+ .uOoOOO 

+.000000 

+.000000 

+. OGOOOO 

+  .OOOOOO 

+. OOOOOO 

+.000000 

+. 472J01 

+  . 36  84  32 

+.71 4bo6 

+.851631 

+.999370 

+.851631 

+.714556 

+.3654 32 

-.265847 

+.000000 

+.OjuOjj 

+. OOOOOO 

+.000000 

+. OOOOOO 

+. OOOOOO 

+. OOOOOO 

+. Oooooo 

+  .  0  0  0  0  u  0 

+  .00  uooo 

+.uOOOGO 

+. OOOOOO 

+ .oooooo 

+. oooooo 

+ .OOCOOJ 

-.41 U003 

- .46  7056 

— . 1 696u 1 

+. 1012G7 

+1 .97931 0 

+.101207 

-. 16960  1 

-.467056 

+.252423 

+.000000 

+  .  0  J  0  0  0  u 

+.000000 

+. OOOOOO 

+. oooooo 

+.000000 

+.000000 

+. OOOOOO 

+.000000 

+  *  0  00  000 

+.000000 

+. OOOOCO 

+ .oooooo 

+. OOOOOO 

+.000000 

+.365432 

+.098474 

+.365432 

+.472001 

+.851631 

+.999370 

+.851631 

+ .47200 1 

-.2471 16 

+.000000 

♦  •  0  u  J  0  J  0 

+.UCGU0G 

+. OOOOOO 

+. OOOOOO 

+. oooooo 

+.000000 

+. OOOoOO 

+  . 0  0 0 0  0 0 

+  . 0  uO uJv 

+. OOOoOO 

+.000000 

+  .OOOOOO 

+.000000 

+. OOOOOO 

- . 467056 

-  .=>2u633 

- . 467056 

-.410003 

+  .  1  J  1  £07 

+1 .979310 

+.101207 

-.41 0003 

+.267421 

+. OOOOOO 

+.OudUUV 

+ . ouoOG J 

+ • OuOOuO 

+. UOOOOO 

+. OOOOOO 

+ • OOOOOO 

+.000000 

+. oOOOOO 

+  .  0  -/O 00 O 

+ . uOouCG 

+. OuOOOO 

+  .OOOOOO 

+ .OOOOOO 

+ • OOOOCO 

+.71 4556 

+  • 3654  32 

+  .4  7200 1 

+  •  36o>432 

+.71 4556 

+.851631 

+.999370 

+.851631 

-.396862 

+ .uOUOUQ 

+ . 0 0 uOuu 

+  .  u^OCiOO 

+ . OOOOOO 

+ .OOOOOO 

+.OOOuOO 

+.000000 

+ . OOOOOO 

+ .OOoOOO 

+  ..  0  0 0  G  j  0 

+ . OOCuOG 

+.000000 

+  .  U  0  0  0  0  0 

+.OCOOOO 

+ .OOOOOO 

-. 169601 

- .467056 

-.4 l 0003 

-.467056 

-. 1 69601 

+.101207 

+1 .97931 0 

+. 101207 

+.018899 

+ .OGoOOO 

+.0000^0 

+ . OOOOOO 

+. JoOOOO 

+  . OOOOOO 

+.COOOOO 

+  .OOOOOO 

+. OuOOOO 

+. OOOOOO 

+  .  0  u  0  00  0 

+ . oooooo 

+. OOOOOO 

+  .oooooo 

+. OOOOOO 

+ .OOOOOO 

+.851631 

+.472001 

+.305432 

+.098474 

+. 365432 

+.472001 

+.851631 

+.999370 

-.329898 

+. OOOOOO 

+»09JvJj 

+.000000 

+. OOOOoO 

+ .003000 

+. OOOOuO 

+  .COOOOC 

+.000000 

+. OOoOOO 

+ .  o  ouoo  «j 

+  .'-'00000 

+. OuOOOO 

+  .uOOGOU 

+. OOOOOO 

+  . OOOOOO 

+. 10 1207 

-.4 1 0003 

- .467056 

-.52^633 

- .467056 

-.4 1 0003 

+  .  101207 

+1 .979310 

-.292726 

+.OCuOGO 

+  «  v;  J  u  9  J 

+  •  > .  u  u  C  u 

+• OOOooO 

+ .uOOOOO 

+. OOOOOO 

+ .OOOOOO 

+. OOOoOO 

+  .uOoO  jO 

+  .  0  J1; 

+  .J'.uuO  J 

+ . uooOOO 

+ .OOuOOO 

+. OOOOOO 

+.000000 

Fig.  20.  Inputs  and  outputs  of  the  programs  in  appendices 
A,  B,  and  C . 
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C  INPUT  CARDS  FOR  PROGRAM  in  FIG  16 

a 


+.999370 

+.851631 

+.71 4556 

+.365432 

+.472001 

+.365432 

+.714556 

+.851631 

-.304243 

+  .000000 

+.000000 

+.000000 

+.000000 

+.000000 

+. OOOOOO 

+.003000 

+ 1.97931 0 

+  .  1  0  1207 

-. 169601 

- .467056 

-.41*0003 

-.467056 

-. 169601 

+.101207 

-.328366 

+  .000000 

+. oooooo 

+.000000 

+.000000 

+. OOOOOO 

+.000000 

+.000000 

+.851631 

+.999370 

+.851631 

+.472001 

+.365432 

+.090474 

+.365432 

+.472001 

-.329698 

+ .000000 

+.000000 

+. OOOOOO 

+.000000 

+.000000 

+.000000 

+.000000 

+. 101207 

+1 .979310 

+. 101207 

-.410003 

-.467056 

-.520633 

-.467056 

-.410003 

-.292726 

+ .000000 

+. OOOOOO 

+. OOOOOO 

+.000000 

+. OOOOOO 

+.000000 

+.000000 

+.714556 

+.851631 

+.999370 

+.851631 

+.714556 

+.365432 

+.472001 

+.365432 

- .396862 

+  .000000 

+.0000U0 

+.0^0000 

+.000000 

+.000000 

+.000000 

+.000000 

-. 169601 

+  .  101207 

+1 .979310 

+. 1 0 i 207 

-. 169601 

-.467056 

-.410003 

-.467056 

+.018099 

+  .000000 

+. OOOoOO 

+.000000 

+.000000 

+ .OOOOOO 

+.000000 

+.000000 

+.365432 

+.472001 

+.851631 

+.999370 

+. B51631 

+.472001 

+.365432 

+.098474 

-.2471 16 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

-.467056 

-.41 0003 

+. 101207 

+ 1 . 979310 

+. 101207 

-.410003 

-.467056 

-.520633 

+.267421 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+. oooooo 

+.000000 

+.472001 

+.365432 

+.714556 

+.851631 

+.999370 

+.851631 

+.71 4556 

+.365432 

- .265847 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

-.41 0003 

-.467056 

-.1 6960 1 

+  .  101207 

+1.979310 

+.101207 

-.169601 

-.467056 

+.252423 

+  .000000 

+ .OOOOOO 

+.000000 

+.000000 

+. OOOOOO 

+.000000 

+.000000 

+.365432 

+.098474 

+.365432 

+.47200) 

+.851631 

+.999370 

+.051631 

+.472001 

-.2471  16 

+.000000 

+. OOOOOO 

+. OOOOOO 

+.000000 

+.000000 

+.000000 

+.000000 

-.467056 

•  -.520633 

— . 467056 

-.410003 

+. 1 01207 

+1 .979310 

+  .  101207 

-.410003 

+.267421 

+  .000000 

+. OOOOOO 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+.71 4556 

+.365432 

+.472001 

+.365432 

+.714556 

+.851631 

+.999370 

+.851631 

-.396862 

+.000000 

+.000000 

+. OOOOOO 

+.000000 

+.000000 

+.000000 

+. OOOOOO 

169601 

-.467056 

-.410003 

-.467056 

-. 169601 

+  .  101207 

+1.979310 

+.101207 

+.010899 

+.000000 

+.000000 

+. OOOOOO 

+.000000 

+.000000 

+.000000 

+.000000 

-r. 851631 

+  .47200 1 

+.365432 

+ .098474 

+.365432 

+.472001 

+.851631 

+.999370 

- .329898 

+  .000000 

+. OoOOOO 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+. 101207 

-.410003 

-.467056 

-.520633 

-.467056 

-.410003 

+.101207 

+1.979310 

-.292726 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

+.000000 

C  OUTPUT  CARDS  FROM  PROGRAM  JN  FIG  16 


+  8 

-.161631 

-.018926 

+  7 

-.038466 

+.028586 

+  6 

+.005762 

+.019949 

+  5 

-.000098 

+.000905 

+  4 

+.005762 

+.019949 

+  3 

-.038466 

+.028586 

+  2 

-.161631 

-.010926 

+  1 

-.133808 

-.062170 

I  nput 

CARDS  FOR 

PROGRAM  IN  FIG  18 

9 

0.0 

1  0.0 

36.0 

o 

• 

o 

90.0 

0.500 

0.500 

0.625  0.750 

0.750 

0.750 

0.625 

0.500 

0.0 

0.0 

O.UUO 

0.  1  25 

0.125  0.125 

0.000 

-0.125 

-0.125 

-0.125 

C  .0 

0.0 

9 

-1.0 

0.0 

+  8 

-.161631 

-.018926 

+  7 

-.038466 

+. 028586 

+  6 

+.005762 

+.019949 

+  5 

-.000098 

+. 00u905 

+  4 

+  .005762 

+.019949 

+  3 

-.038466 

+.028586 

+  2 

-.161631 

-.018926 

1  1 

-.1 33808 

-.0621 70 

Fig.  20.  Inputs  and  outputs  of  the  programs  in  appendices 
A,  B,  and  C.  7 
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C  OUTPUT  CAWDS  FWOM  PROGWAM  IN  FIG  18 


+  9 

-  1 .OOOOuOOO 

+ • OoOO-jCOO 

+ 1 . OOOuOOOo 

+ 1 .OOOOOOOO 

+  8 

16163100 

-.01892603 

+.02646277 

+. 16273528 

+  7 

-.03846600 

+ . 026586u0 

+.00229679 

+  .  04  792466 

+  6 

r. 00676200 

+  .0  1  99490.0 

+.000431 16 

+  .02076447 

+  5 

-.00009800 

+ . O009360O 

+.00000082 

+.00091029 

+  4 

+.00576200 

+.01994900 

+.000431 16 

+.02076447 

+  3 

-.03646600 

+.02858630 

+.00229679 

+.04792486 

+  2 

-. 16163100 

-.01892600 

+.02648277 

+. 1 6273528 

+  1 

-. 13380800 

-.06217000 

+.02176969 

+. 14754555 

+  1 

+.00000000 

, -.407301 16 

♦  .  1 0256843 

+. 1 7641451 

+.4200 1 727 

♦  2 

+10.00000000 

-.40757072 

+.07192805 

+. 17128753 

+.41386897 

+  3 

+20. 00000000 

-.4171 7391 

-.01696248 

+.17432180 

+.41751862 

+  4 

+30. OOOOUOOO 

-.45908866 

-.15166895 

+.23376587 

+.48349340 

+  5 

+  4  c • oooooooo 

-.55959657 

-.30378805 

+ .40543550 

+.6367381  7 

+  6 

+bo. uooaoouo 

-.73096286 

-.42788042 

+.71738835 

+.84698780 

+  7 

+60. 00003000 

-.9554 1 772 

-.47370861 

+ 1 . 13722310 

+  1 . 0664  0660 

+  8 

+70.00000000 

-1.1 8309290 

-.40961800 

+ 1 .56749570 

+ 1 .25199660 

+  9 

+60. OOOOoOuO 

-  1 .35024300 

-.24259704 

+ 1 .88200950 

+1 .37183350 

+  10 

+90. OOOOCOOC 

-1  .4  0873-720 

-.01 939185 

+ 1 .98491650 

+1.40887060 

+  1  1 

+100.03000000 

-  1 .34630170 

+.19377559 

+ 1 .85546650 

+ 1 .36215500 

+  12 

+1 1 0. oooooooo 

-1.1 9648003 

+.34034474 

+ 1 .552 18880 

+ 1.24586860 

+  13 

+ 1 20.00000000 

-1.011 33640 

+.39454 1 28 

+1.1 78464 1 0 

+ 1 .08557080 

+  14 

+130.C000O00U 

-.83707761 

+.36484096 

+.83380786 

+.91313030 

+  1  5 

+140.00000000 

-.73644604 

+.28205393 

+.57862043 

+.76067103 

♦  16 

+ 150. OOOOOOOO 

-.62652430 

+. 18230454 

+.42576765 

+.65250B73 

+  17 

+160. OOOOuOOO 

-.5874 1117 

+.09461487 

+.35400386 

+.59498223 

+  18 

+170. uOOOJOOO 

-.5729231 7 

+.03666559 

+.32958533 

+.57409522 

+  19 

+ 1 90. OOOOoOuO 

-.56976057 

+.01562214 

+.32490341 

+.57000299 

+  20 

+190. JOOOOOJO 

-.57292313 

+.03666538 

+.32958526 

+.57409516 

+  21 

+2uO. 0000 JOOO 

-.5874 1 1 06 

+.09461452 

+.35400366 

+.59498206 

+22 

+210. OOOCOOOO 

-.62652402 

+.18230404 

+.4257671 1 

+.65250832 

+  23 

+220.00000000 

-.70644550 

+.28205340  ' 

+.57861936 

+.76067033 

+  24 

+230.00000000 

-.83707676 

+.36484062 

+.83380618 

+.91312988 

+  25 

+240.0000 OOJO 

-1.01133550 

+.39454 1 30 

+1.1 7846230 

+ 1 .08557000 

♦  26 

+250. OOOOOOOo 

-l . 19647910 

+.34034526 

+1 .55218710 

+ 1 .24586800 

+  27 

+260. 0000 JOJO 

-1  .34830  1  10 

+.19377658 

+1 .85546530 

+1 .36215460 

+  28 

+270. OOOCOOOO 

-1.40873720 

-.01939065 

♦1.98491650 

+1.40887060 

+  29 

+280.00000000 

-1 .35024360 

-.24259598 

+1 .88201060 

+ 1 .37186390 

+  30 

+2yu • ouco  ooo o 

-  1  .  1 B3O9400 

-.40961738 

+ 1 .56749780 

+1 .25199750 

+31 

+300. 0000 JOOO 

-.95541896 

-.47370878 

+1.1 3722540 

+ 1 .06640770 

+  32 

+3 1 0. OOOOUOOO 

-.73096395 

-.42788091 

+.71739037 

+.84698900 

+  33 

+320.00000000 

-.5595972 7 

-.30378879 

+ .40543673 

+.63673913 

+  34 

+330. OOOOOOOO 

-.45908900 

-.15166970 

+.23376641 

+.48349396 

♦  35 

+340.00000000 

-.41717405 

-.01696307 

+.17432194 

+.41751878 

+  36 

+350.00000000 

-.40757074 

+.07192772 

+.17128751 

+.41 386895 

Fig.  20.  Inputs  and  outputs  of  the  programs  in  appendices 
A,  B,  and  C . 


41 


APPENDIX  D 


Figure  21  is  a  SCATRAN[4]  program  for  use  on  the  IBM  7090 
computer  at  the  Ohio  State  University.  By  combining  the  programs 
in  Appendix  -B  into  a  single  SCATRAN  program  and  utilizing  the  greater 
speed  of  the  IBM  7090,  this  program  achieves  the  same  results  in 
about  l/50th  the  time  it  takes  the  IBM  1620,  In  addition,  the  IBM  7090 
provides  more  storage  capacity  which  is  necessary  for  the  solution  of 
arrays  with  large  numbers  of  wires. 
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THIELE. 
SA 
F  SI 


F  S2 


S4 

F  S5 
F  SlO 

S12 

SI  4 

SI  6 


G.  J08N  FGJ176  01/13/65 

WRITE  OUTPUT  .SI-  0026 

<///78H  SCATTERING  BY  ARBITRARY  ARRAY  OF  THIN  WIRES  WITH  A  0029 

LINE  SOURCE  AT  THE  ORIGIN///)  -  0030 

I JK=62—  '  0031 

DIMENSION  <U(61).V(61>)-  0032 

DIMENSION  ( A (3782. 1 JK > . B < 3782  * I JK ) )-  0033 

(7F10.6)  -  0034 

READ  INPUT  .S3 .(III)-  0033 

DO  THROUGH  ( S260 ) . JJJ= 1 . 1 . J J J. LE • I  I  I -  0036 

READ  INPUT  .S2. <FN.BKJ.BKY.AAA.PHI )-  0037 

M=FN-  0036 

READ  INPUT  .S2. < <U< I  )  .  I  *1 . 1 . I .LE.N)  1-  0039 

READ  INPUT  .  S2> <  (V (  1  > . I *1  *  1 . I .LE.N))-  0040 

DO  THROUGH  t S4 ) .  I  *  1 . 1 . I .LE.N-  •  004  1 

U< I )=6.2831853*U< 1 )-  0042 

V< l  >=6.2831853*V<  I  1-  0043 

NN=N+1—  0044 

WRITE  OUTPUT  .510-  ’  0045 

(2I5.F15.8)  -  0046 

I///27H  INTERWIRE  SPACING  KR(!.J)///|  •  -  0047 

DO  THROUGH  C S50 ) . I = 1 . 1 , I .LE.N-  0048 

DO  THROUGH  <S30).J*1 .l.J.LE.N-  0049 

TRANSFER  TO  < Si 4 . S I  2 . S 1 4 )  PROVIDED  (I-J1-  0050 

A(I«J)=3KJ-  0051 

B(I«J)*-BKY-  0052 

TRANSFER  TO  (S301-  0053 

R  =  SORT . <  <UC I  )-U< J) >*<U< 1  >-U  < J  > )  +  < V I  1  )  —V ( J ) ) * < V ( I  >-V ( J ) ) >-  0054 

WRITE  OUTPUT  .SS.U.j.Rl-  0055 

TRANSFER  TO  < SI 6»S 1 6 . S 1 8 )  PROVIDED  <R-4. )-  0056 

W=R/2 • —  -  0057 

W 1 =W*W-  0058 

W2=W1*W1-  0059 

W3=W2*W1-  0060 

W4«W2*W2-  '  0061 

W5«W4*W1-  0062 

W6=W3*W3-  0063 

W7=W3*W4-  0064 

B8»l .-Wl  +  .23*W2-.0277778*W3+.00 I  7361 *W4-  0065 

BB »88-. 6944.X. -4 *W5+. 1 929 . X . -5*W6- . 394 . X . -7*W7-  0066 

S*W1-.375*W2+.0509259*W3-.0036169*W4-  0067 

S*S+« 15856.X. -3*W5-. 4726. X.-5*Wf +. 1021 .X.-6*W7-  0068 

Y  =  *636619774  <  S+BB*  (  .  5772 1566+LN.  (Win-  0069 

TRANSFER  TO  <S20)-  1  0070 


Fig.  21.  SCA.TRAN  program  for  chaining  the  far-field  pattern 
of  an  arbitrary  array  with  a  line  source  at  the  origin. 


SIS  Z=l./R-  0071 

X=Z*Z-  0072 

P« ( 1 .-.07031*X+. 1 1215*X*X 1*. 5641 896-  0073 

0  =  < . 125-.07324*X 1  #Z* . 564 l 896-  0074 

BB*<  (P-Q)*C0S. <Rl  +  <PtQ)*SIN.(R ) 1/SORT. (R1-  0075 

Y»  (  (P-Q  1 *SIN.  IRl-lP+O 1*C0S.  ( R  )  l/SQRT.  (R  )-  0076 

S20  All  *  J  >  *88—  0077 

8<1«U1=-Y-  0078 

S30  CONTINUE  -  0079 

R«SQRT.  lull  )*U(  I  1+Vl  IXVIIII*  0080 

TRANSFER  TO  < S36 • S36 • S3B )  PROVIDED  (R-4.1-  0081 

S36  M=R/2.~  0082 

M1=W*M-  0083 

M2 = W 1  *  W 1  -  0084 

M3«W2*W1-  0085 

M4=W2*W2-  0086 

«5=M4*W1-  0087 

M6»W3*W3-  0088 

M7=W3*M4-  0089 

88=1  .-Ml  +.25*M2-.0277778*M3+.001  7361*M4—  0090 

88=00-. 6944.X. -4*MS+. 1  929  •  X  .  -  a*  M  6- •  394  .  X  .  -  7*»7-  0091 

S=W  1  -  .375*1*2  + »  0509259  *M3-»  0036  169«M4-  0092 

S=S+. 15856. X.-3»M5-.4726.X.-5*M6+. 1021 .X.-6*M7-  0093 

.  Y=.e3661977#<S+8B*< . 5772 1 566+LN.  tM >  1  >-  0094 

TRANSFER  TO  (S491-  0095 

S38  Z=l./R-  0096 

X=Z*Z-  0097 

P=( 1 .-.07031 *X+. 1 1215*X*X )«. 5641896-  0098 

0= ( . 125-.07324*X)*Z*.5641896-  0099 

88= ( <P-Ul*COS.<Rl+(P+Ql*SIN.<Rl l/SQRT. (Rl-  0100 

Y=( <P-Q|*SIN. <R )-<P+Ol*COS. <Rl l/SQRT. (R  1-  0101 

S49  B ( I tNN 1 =— Y—  0102 

S30  A<I.NN1=8B-  ~~  0103 

MR1TE  OUTPUT  . S2 . I < I  A ( 1 . J 1 ♦ J* I ♦ 1 . J.LE.NN 1 . 1 » 1 « 1 • I *LE .N I  1 -  0104 

MR 1TE  OUTPUT  • S2 . ( ( ( 8 ( I . J 1 . J« I . I . J.LE.NN 1 . I » 1  *  1 . 1 .LE.N ) )-  0105 

MR1TE  OUTPUT  ,S60-  0106 

F  S60  I///46H  REAL  AND  IMAGINARY  PARTS  OF  CURRENT  ON  MIRE  !///>  0107 

0108 

F  S80  ( 1 5 • 2F 1 5 • 8 )  -  0109 

00  THROUGH  < S l 1 8 1 .L *  1 . 1 .L.LE .N-  0110 

LLL  =  L  -  1  - 

DO  THROUGH  1  SI  18 1*1 =L. 1*1 .LE.N—  0112 

1 1*1+1—  0111 

TRANSFER  TO  (S 1 04 *S1 04 * S 1 03 )  PROVIDED  (LLL 1—  0114 


Fig.  21.  SCATRAN  program  for  obtaining  the  far-field  pattern 
of  an  arbitrary  array  with  a  line  source  at  the  origin. 
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3103  00  THROUGH  ( 5 1 1 7 > . K« 1 . 1 . K .LE .ILL-  0115 

A(  I  .L ) =A ( I 4L ) -A  <  1  .X  >*A (K.L  >+3 ( I  .K )#B (K.L>-  0116 

B<  I .L ) =B  (  l  . L)-A(l. K  >*B  (K  .L  ) -A < K . L  > *B  <  l . K  >-  0117 

A  <L« !  I  )  =  A(L.  1 1  >-A <L»K  >*A (K> I  1  )+B(L.K >*B(K.  I  I  >-  01  IB 

S 1 1 7  B (L. 1  1  >=B (L.  I  1  )-A<L.K >*B (K. 1  I  )-8<L.K )*A<K.  I  I  )-  0119 

S 1 04  R=A (L.L )*A< L.L )+B (L.L )*B (L.L )-  0120 

S=A(L.I11-  0121 

T  =B ( L  « I  1  )-  0122 

AIL* I  I  1= < A(L<L )*S+8 (L.L  )*T  l/R-  0123 

sue  aiLtin«(A(Ltu*T-B(Ltu*si/R-  0124 

WRITE  OUTPUT  <  580  .  ( N «  A INiNN )  •  B ( N« NN ) ) —  0125 

DO  THROUGH  ( S 123 ) • L=2 < 1 «L .LE 0126 
I=NN-L-  0127 

11=1+1-  0128 

DO  THROUGH  ( S 1 22 ) . K= I  I « 1 . K.LE.N-  0129 

A (  I  <NN ) =A  < 1 . NN )-A ( I  .K  )*A IKiNN >  +  8f  I *K )#8(K. NN )-  0130 

5 1 22  B( 1 .NN)=B( I «NN)-A( 1 . K ) +B < K . NN )-8 ( 1 .K )*A 1K» NN ) -  0131 

5 1 23  WRITE  OUTPUT  . S8Q « ( I . A ( I . NN ) «B ( 1 . NN 1 >-  '  0132 

WRITE  OUTPUT  .S200-  0133 

F  5200  1///18H  FAR  FIELD  PATTERN///)  -  0134 

READ  INPUT  . S2« tP.DP.FM )-  0135 

M=FM-  0136 

A ( 0 . NN 1  *  —  1 —  0137 

B ( 0 . NN ) * . 0—  0138 

U<0)=.0-  0139 

V(0)=.0-  0140 

DO  THROUGH  (S240). 1«1 .1 . l.LE.M-  0141 

PR= . 0 1 745329+P-  0142 

C  =  C0S. (PR  1—  0143 

S*S 1 N. <  PR ) -  0144 

ER=0.-  0145 

El =0.-  0146 

DO  THROUGH  ( S230 ) • U=  0 « 1 .J.LE.N—  .  0147 

2=U< J)*C+V< J)*S-  0148 

CZ=COS.<Z>-  0149 

SZ  =  S 1 N. <  Z I —  0150 

ER=ER+A ( J.NN ) *CZ— 8 ( U .NN ) *SZ—  0151 

S230  E1=EI+B( J.NN)*CZ+AC J.NN)*SZ-  0152 

ES=EW*ER+E 1 *E I  -  0153 

E=S0RT.<ES>-  0154 

WRITE  OL'TPUT  .5250.  tP.ER.EI  .ES.E  )-  0155 

S24G  P=P+DP-  0156 

F  S250  ( 5F 14.8)  —  0157 

S260  WRITE  OUTPUT  ,51-  0158 

END  PROGRAM  ( SA )—  0159 


Fig.  21.  SCA TRAN  program  for  obtaining  the  far-field  pattern 
of  an  arbitrary  array  with  a  line  source  at  the  origin. 
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